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1. Introduction 

Heavy flavor physics invariably plays an important role in the lattice community. As with the 
light quark sector, quantities involving heavy quarks can be studied to both test the methodology 
used in lattice calculations as well as testing the physics of the Standard Model. Depending on 
the situation, if there is a discrepancy between experiment and a lattice result, one must be able 
to determine the source of this discrepancy, and whether or not it actually exists. Thus, whether 
one is testing the lattice techniques or the Standard Model, reliable comparison with experimental 
measurements is crucial. 

Reliable comparison with experiment takes two forms. First, one must have good quantitative 
control over all systematic uncertainties which enter into a calculation, giving credence to the 
results. However, there is another criterion which is implied but rarely explicitly stated, even when 
stressed as important in other arenas, and that is agreement among different calculations of the 
same quantity using different formulations. This is an area where the light quark sector (especially 
in kaon physics, see Vittorio Lubicz's review at this conference [[I]]) has a tremendous advantage 
over the heavy quark sector. This becomes important in cases such as the "/d s puzzle," where a 
significant discrepancy between experiment and the lattice calculation has shown up, and for two 
years has not yet gone away. This discrepancy exists between the experimental average and a single 
lattice calculation, and has not be reproduced by another. 

For this review, I divide my discussion into three parts. 1 I begin with a couple of examples 
of spectroscopy calculations which simultaneously allow lattice calculations to be tested as well 
as predicting experimental measurements. Specifically in Sec. ||, I discuss one meson example 
(that of the B* mass) and the single-bottom baryon spectrum, where there has been some interest 
during the last year when comparing lattice and experimental results. Then I will discuss cases 
of calculations which are more focussed on testing the Standard Model. In Sec. & I will discuss 
current results on B-mixing and decays, which is at a stronger point than most as it has a larger 
number of groups which allow for decent cross-checks. In Sec. ||, I will discuss D decays, mostly 
sticking to a discussion of the possible discrepancy between the lattice calculation and experimental 
measurement of fo s . In the final part, I discuss newer ideas and what lies in the future for heavy 
quark physics. In Sec. ||, I discuss new ideas for extracting semileptonic form factors for D decays, 
and in Sec. |6[ I present some initial ideas for attempting extractions of nonleptonic B decay matrix 
elements from a lattice simulation. I finish with some conclusions in Sec. 0. 

2. Spectroscopy 

Lattice studies of the particle spectrum can play two roles in understanding Nature. First, 
it is an invaluable tool to test lattice techniques by calculating particle masses that are known 
experimentally to see that the methodology used is sound. At the same time, it can be use to 
predict masses for states that have not yet been seen. One nice feature in spectrum calculations 
that involve heavy flavors is that often lattice calculations have been done at roughly the same time 

'Obviously I cannot provide an exhaustive review of what has been reported in the field during this past year. I will 
mostly focus on a few select results, and only on those with at least two flavors of dynamical light quarks. I apologize to 
all those whose calculations I didn't have time to discuss. 
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Figure 1: Plot of niB* versus the lattice spacing squared, showing the extrapolation to the continuum limit. 
On the same plot is the lattice determination of ms c from [g, ||| compared with the PDG Q. Plot courtesy 
of E. Gregory, from [||] . 



as the experimental measurements, so there are more cases of predictions, as opposed to the usual 
postdictions, of the particle spectrum. 

An example of this is the measurements of the B c system by the HPQCD collaboration. In 
[Q, 0, using NRQCD b quarks and HISQ c quarks in the valence sector and the MILC 2+1-flavor 
AsqTad configurations (three lattice spacings), HPQCD predicted the value of the B c mass, shortly 
afterwards to be confirmed by experiment [Bp. This year, they have determined a prediction for the 
vector meson, the B* using the same techniques. In this case, they actually measure the ratio 

mm — mn 

^ , (2.1) 

m B * - m Bs 

and use the experimental value (mg» — m#J exp = 49.0(1.5). The results are shown in Fig. [I], and 
they obtain mg» — mg c = 0.059(6) GeV for the mass splitting [Q]. 

On the baryon side, there have been many lattice results of baryons with one or more bottom 
quarks [[7], ||, ^, [K], 11]. Of these, there has been some interest in the last year in the single-bottom 



sector, due to an apparent discrepancy between lattice results and experiment. 

A summary of some results is shown in Fig. |2|, where I compare the lattice results 2 to the 
experimental results from the DO collaboration [ |I2"| , [T3[ ] and the CDF collaboration JT4] , [T5| , [16| ]. 
The puzzle as of last year was the agreement between the results shown in red (a calculation using 
NRQCD heavy quarks and Clover light quarks by Lewis and Woloshyn, [[7|]) and DO for all masses 
except for the Q.},, where there was an unexplained difference between the DO measurement and 



2 Note that only results which use 2+1-flavors of sea quarks are shown, and all results to date have only been 
presented at a single lattice spacing. 
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Figure 2: Singly-bottom baryon spectrum from various lattice calculations and experiment. The data in 
red are using NRQCD heavy quarks [Q], in blue use static heavy quarks Jl0|], and the green data are from 
experiment [|(| [j~2|, [lj, [l5| [j~3| Efl]. The measurement of the mass of the Q./, differs between DO (black) and 
CDF (green). Data from Refs. [[j (static heavy quarks) and [ 11 1 (NRQCD quarks) are not shown for clarity, 
but is consistent with the other lattice calculations. The vertical scale is in MeV. 



theory. Later, various results using static heavy quarks |g[ [lO|] and NRQCD quarks [11] confirmed 
the earlier lattice results, also shown in Fig. ||. 3 

The results from CDF published in May fll6[ ] confirm the lattice picture, and thus the dis- 
crepancy has shifted from a lattice-experiment discrepancy to one between differing experiments, 
which I will not comment on further, as this is not within the scope of this review. What is impor- 
tant to note here is that a single lattice calculation which disagrees with experiment is not evidence 
enough for new physics, nor can we use this to claim either the experiment or the lattice calculation 
is flawed in some way. Rather, only when there are multiple lattice calculations (and experimental 
measurements) can we claim (dis)agreement between the the different sides. 



3. 5-physics 

Moving on to calculations of matrix elements, I begin with the &-quark sector. Here, there 
have been calculations for the B and B s meson decay constants as well as for the neutral B s and 
Bd mixing parameters. Often these calculations are performed in concert with each other, because 
if one is interested in the fi-parameters themselves (ie, separated from the decay constants), one 
needs the decay constants. 

As is well known, the decay constant fg governs the leptonic decay rate of a B-meson to a 
lepton-antineutrino pair. The B-parameter is defined by (B q \OLL\B q ) = |/J ms Bb ,defined in 
this manner so that Bb ci = 1 in the vacuum saturation approximation, and Oil is the left-left weak 
operator governing this transition. The B% q parameter enters explicitly into the expression for the 

3 The lattice results are consistent with theoretical understanding coming from Heavy Quark Effective Theory and 
various quark model pictures. 
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2+1 flavor 



Group 



a 
(fm) 



m 



JT.mm 



Q 



f B ,f B , B B Refs 



(MeV) 



Fermilab/MILC 
HPQCD I 
HPQCD II 

RBC/UKQCD 



0.09,0.12 230 

0.09,0.12 260 

0.09,0.12,0.15 320 

0.11 400 



Asqtad Fermilab / 

Asqtad NRQCD / 

HISQ NRQCD 

DW Static / 



/ PIPD 

/ 



[19] 



[20] 



[21] 



2 flavor 


Group 


a 
(fm) 


(MeV) 


q 




Refs 


ETMC 


0.05, 0.065, 


300 


TM 


Static/TM / 


[ 


22 


] 



Burch et al 



0.085,0.10 
0.11,0.16 



350 



CI 



Static 



Table 1: A summary of the various calculations of the B,B S meson decay constants and mixing parameters. 
q denotes the light quark action (where DW = Domain-wall, TM = Twisted Mass Wilson, and CI = Chirally 
Invariant), and Q denotes the heavy quark action used. A checkmark denotes that the quantity (or at least the 
ratio of quantities) have been reported before or during this conference. 



oscillation frequency of the B q meson, given by 

AMq = ~6? I VtqV t b | 2 rj^So (x t )m Bq fl q B Bcj , (3.1) 

where the left-hand side is determined experimentally, and only the decay constant and B Bq on the 
right-hand side are needed nonperturbatively. What is generally of most interest phenomenologi- 



cally is the CKM matrix elements in Eq. ( |3.1[ ), specifically the ratio 

\Vtd\ 




\ v ts\ fB^/BZv km* m B 



(3.2) 



The nice feature of this ratio is that only the quantities that are directly measured either on the 
lattice or experimentally are required to obtain this ratio. The quantity from the lattice is the ratio 
of the mixing parameters and decay constants, and is denoted by £ = (f Bs ^B Bs ) / (f Bd \jB Bi ) . 
Additionally, as with any ratio, many systematic errors will cancel, and as such a more precise 
determination of £, can be made compared with any of the quantities individually. 

In Table |l] I summarize the various collaborations who have calculated either the decay con- 
stants or the mixing parameters, or both (or in some cases just ratios of them). I list the lattice spac- 
ing^) used, the minimum pion mass, as well as the formulations for the light and heavy quark ac- 
tions. All simulations use the same valence and sea light quarks, except for that labelled "HPQCD 
II," which uses HISQ valence light quarks and AsqTad sea quarks, and this is the main difference 
between the two HPQCD calculations listed. Additionally, while there are preliminary results for 
the HPQCD II calculation, no numbers have yet been presented resulting from this calculation. 

The Fermilab/MILC collaboration has not updated their results for the B decay constants this 
year, however they have for their determination of B B , or more precisely, ^ . As discussed in the 
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2+1 flavor 


Group 


fs (MeV) / Bl (MeV) 


/bJ/b 




Fermilab/MILC 


195(11)* 243(11)* 


1.245(43)* 


1.205(50) 


HPQCD I 


190(13) 231(15) 


1.226(26) 


1.258(33) 


RBC/UKQCD (APE) 




1.20(14) 


1.187(112) 


RBC/UKQCD (HYP2) 




1.19(16) 


1.18(19) 


2 flavor 


Group 


fs /b s 


fsJfB 




ETMC 


203(17) 247(16) 


1.22(6) 




Burch et al (0.16 fm) 




1.108(29) 




Burch et al (0.11 fm) 




1.089(41) 





Table 2: A summary of results from the various collaborations for the B decay constants and mixing param- 
eter. The asterisk on the Fermilab/MILC numbers denote that these are the results have not been updated 
since Lattice 2008. For the Burch et al calculation, no continuum limit was taken for the 2-flavor results, so 
I list them separately. 



parallel session [18], the group has checked their results for the perturbative renormalization, and 
as such has a more well-defined determination of the systematic errors. For the ratio E, , this is a 
small percentage (~ 0.2%), and not the dominant error. Most of the uncertainties tend to cancel 
in the ratio, and the only two errors that remain dominant are the light quark discretization and 
chiral fits (using staggered light quarks, they implement staggered XPT fits [^3|, |4|, 25]), and the 
statistical, which are 2.8% and 3.1%, respectively. One can see in Table §that these two dominate 
the total error in £ , which is roughly 4.3%. 

The RBC/UKQCD calculation, being an initial study, has fewer data than the Fermilab/MILC 
group and only a single lattice spacing. However, while the formulations used are different, both 
the central values an sources of the largest systematic errors are similar. They use static quarks for 
the b quark, using two different smearings (APE and HYP2) to control the heavy quark discretiza- 
tion errors, and domain wall light quarks for both the valence and sea. After performing the chiral 
fits, their errors too are dominated by statistics, discretization effects, and XVT. Of the two smear- 
ings, they find smaller systematic errors coming from APE smearing, of a combined 5.4%+7.7% 
(stat+sys) on £ for the APE smearing compared with 4.0%+15.1% with HYP2 smearing. A similar 
(although not as great) reduction of errors is seen for the ratio of /bJ/b, and both are shown in 
Table g. 

As for the 2-flavor simulations, there exist currently only calculations of the decay constants, 
and not £ . ETMC performed simulations using both static quarks and twisted mass Wilson quarks 
for the b quark on four lattice spacings, controlling the continuum limit quite well. Their chiral 
fits are not using the relevant XPT for tm Wilson quarks, but instead are polynomials in the heavy- 
light mass, rrihq, and 1/m^. As for the calculation of Burch et al, their primary focus was on a 
determination of the excited hadron spectrum involving heavy quarks using a variational approach 
with chirally invariant light quarks and static heavy quarks. This allows one to obtain ratios of 
the decay constants quite easily, although they have not (for the dynamical case) performed a 
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Figure 3: Preliminary results from the HPQCD collaboration for fs s ^/niB s using HISQ valence quarks, 
compared to using Asqtad. Figure courtesy of J. Shigemitsu, and presented in [p0||. 



continuum extrapolation, and as such the systematics are incomplete. Both of these collaborations' 
results are shown in Table |2[ 

One can see in Table || that the results from the different collaborations are quite consistent 
for the various quantities. What is useful here is that there are various caculations of the same 
quantities, a couple of which having similarly sized (and well-controlled) uncertainties. Of course 
these quantities are also of tremendous importance to connect with experiment and especially to 



constrain CKM matrix elements [26] 



In addition, HPQCD has repeated their calculation of B-related quantities replacing the va- 
lence AsqTad quarks with HISQ quarks (still using NRQCD for the b quark), and have presented 



preliminary results at this conference [gOJ. As there are no final results yet, I will only show some 
indicative results for fg s , in Fig. |5[ This plot shows fs^fng, as a function of the valence quark 
mass for two (three) different lattice spacings for the AsqTad (HISQ) calculation, as well as the 
continuum extrapolated AsqTad results. One can see very clearly that difference in discretization 
errors between the AsqTad and HISQ calculations, which again shows how the HISQ action re- 
duces the lattice spacing errors significantly. This is quite important as we have seen that lattice 
discretization errors tend to be one of the more dominant contributions to systematic uncertainties. 
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4. f Ds and \ V CS \ 

In the charm quark sector, there has been an ongoing "puzzle" for the last couple of years 
with regards to the D s decay constant. This puzzle was introduced when the HPQCD collaboration 



published an extremely precise result (~ 1.2%) for fo s Q27Q that disagreed by roughly 3a from 



experimental results (J28J, [29J, |30J]. This disagreement is enhanced by the fact that all other quantities 
calculated using the same methodology [HISQ valence quarks (u,d,s,c) on AsqTad sea quarks] 
agree quite well with the experimental measurements. 

To discuss this possible signal for new physics, I first outline the ingredients. The decay rate 
for D s — > £v is given by 

r (D ,^v) = ^ 1 ^i--f) fl\v a \\ <4.i) 

with the lattice being able to provide fo s , and experimental measurements of the decay rate can 
determine only the product /b s |V C s| by measuring T. Assuming the CKM matrix is unitary, we 
have the relation \V u d\ = \V CS \ + ^(A 4 ), where A 4 ~ 0.003 and \V u d\ is known extremely well, so 
one can use this substitution with the experimental measurement to extract fo s . 

As of last year's conference, the discrepancy remained (see E. Gamiz's review from Lattice 
2008 |H]]). The results were 4 

HPQCD: f Ds = (241±3)MeV 
CLEO + Belle: f Ds = (270.4 ± 7.3 ± 3.7) MeV . 

Here I do not list the value determined by the Fermilab/MILC collaborations [|T7j], which was 
consistent with HPQCD with a factor of three larger total uncertainty (and as such, it was also 
consistent with experiment). 

During the last year, there has been significant progress both on the theoretical side as well as 
experimental. CLEO-c has published new results, which I show in Table |3| In this table, adapted 



from Ref. []30|], I show the most recent results from Belle and CLEO-c which determine the absolute 
branching ratios for D+ — > n + v and Df — > T + v. 5 

One can see that the average of the CLEO-c and Belle numbers has come down significantly, 

to 

f Ds = 26 1.2 ±6.9 MeV. 

The error here is also reduced some, so the discrepancy between this and HPQCD is around 2.5 a 
now, which is reduced but possibly significant. 

HPQCD has performed more simulations, adding two lattice spacings (using the superfine and 
ultrafine MILC lattices, a « 0.06 fm and a « 0.045 fm, respectively). The analysis has not been 
complete, but the new data points are shown in Fig. ||(a), with the previous three lattice spacings 



4 I listed only the averages of the experimental results which measure absolute branching ratios for this decay, which 
is not what is listed in the Particle Data Book. The reason for this is to remove the extra (possibly large) systematic 
coming from the relative branching ratio measurements. 

5 Ref. [ {3(]| l includes results that measure branching ratios relative to D s — » Again, I do not consider these due to 
the large systematic uncertainties that plague these measurements. 
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Table 3: Summary of experimental results for &(D+ — > M + v), &(Df —> T + v), and/ D +, taking into account 
results from absolute branching fractions only. Results have been updated for the new value of the D s lifetime 
of 0.5 ps (This table adapted from Table V of Ref. pfl.) 

Exp. 



Belle [28] 



CLEO-c []29p 
CLEO-c 
CLEO-c @] 
CLEO-c 



Mode 








fot (MeV) 


;U + V 


(6.38 ±0.76 ±0.52) 


10" 


-3 


274±16±12 


T+V 


(5.30 ±0.47 ±0.22) 


10 


-2 


252.5± 11.1 ±5.2 


/I + V 


(5.65 ±0.45 ±0.17) 


10 


-3 


257.3 ±10.3 ±3.9 


T+V 


(6.42±0.81±0.18) 


10 


-2 


278.7± 17.1 ±3.8 




combined all CLEO-c results 




259.5 ±6.6 ±3.1 



Average of CLEO and Belle results above, radiatively corrected 



261.2 ±6.9 



as well as the 2007 result (which only involves a chiral/continuum extrapolation from the heavier 
three lattice spacings, not the new data). One can see that the new points fall right in line with 
the other three, and as such at first glance this will not change the final central value much, only a 
reduction of errors will emerge. 

However, one thing that has not been included is a new determination of the relative scale, 
r\ . Since r\ enters the calculation in a highly non-trivial way, it is difficult to know how a new 
determination will affect the final results. As such, it is hard to predict what the updated value for 
fo s will be with the new HPQCD data and inclusion of the the new r\ value. 

In addition, new preliminary results have been presented at this conference from the Fermi- 



lab/MILC [|32j and #QCD [ p3| ] groups, as well as a 2-flavor result from ETMC, which was pub- 
lished in Ref. [34]. The Fermilab/MILC result is an update of previous results using AsqTad light 
quarks and Fermilab heavy quarks on three lattice spacings. The #QCD calculation uses the Over- 
lap formulation for both the light and charm valence quarks on the RBC/UKQCD Domain- Wall 
Fermion lattices. This calculation is done on two lattice spacings, with only fifty configurations. 
The ETMC calculation uses twisted-mass quarks for charm quark as well as the light quarks and 
they have performed the simulation on three lattice spacings, but only have two dynamical quark 
flavors. 

I show the various 2 and 2+1-flavor results to date in Fig. f|(b). The latest 2-flavor result from 
ETMC [34] agrees with the previous results from several years ago [[35|, 36], with significantly 
improved errors. The two Fermilab/MILC numbers are from those reported at Lattice 2008 [17] 
and a newer number presented at this conference [32]. One thing to note is that except for the 
#QCD calculation [33], all numbers shown on the plot have included both statistical and systematic 
uncertainties (added in quadrature), while #QCD currently only quotes a statistical error. Another 
interesting point is that the new Fermilab/MILC number has shifted upwards, slightly outside (but 
not significantly) the errors of the HPQCD calculation. 

Before making any conclusions, however, I would like to point out that CKM unitarity plays 
a vital role in the experimental determination, and the question remains as to whether or not the 
relation \V u d\ ~ \V CS \ actually holds, and the lattice can play a part in understanding this. This 
involves calculating the semileptonic form factors for D s — > Klv to determine |V CS |, which has 
been done for 2+1 -flavors by the Fermilab/MILC group (initially presented in Ref. [|37|]), with the 
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Figure 4: (a)Updated HPQCD results including the superfine and ultrafine MILC lattices. The black point 
is the extrapolation including only the three original data points (filed circles), (b) All results with either 2 
(red) or 2+1 (blue) flavors of dynamical quarks for fo s , compared with both the PDG value as well as the 
most recent CLEO-c result. The larger symbols (for ETMC, Fermilab/MILC '09, and ^QCD) denote lattice 
results presented within the last year. The dashed error bars on the ^QCD result indicate that these errors 
include only statistical uncertainties. 



most recent determination presented last year being 

|V C ,| = 1.015±0.015±0.106 . (4.2) 
This is of course consistent with unitarity, but with an 1 1% error, this is not a stringent constraint. 



New/updated calculations by ETMC [38], Fermilab/MILC [32], and the Regensburg group [39] 
have been discussed at this conference, but no results for the CKM matrix element has yet been 
presented. 

So the question remains, is there really an fj) s puzzle? There is no conclusive evidence that 
there is a discrepancy between the lattice calculations and experiment as yet. This emphasizes 
the need for other determinations of fo s to match the precision quoted by HPQCD. As it stands 
now, it may be that the HPQCD is an outlier, and this discrepancy may slowly disappear (as these 
"new physics" indications often do). However, with more calculations improving, this could be an 
indication of new physics in the charm quark sector. 



5. New Methods for semileptonic decays 

As discussed at the end of the previous section, there is an increasing need for precision in 
calculations of semileptonic decays of D — > Plv, where D is either a D + or D s meson, and P can 
be % or K. The difficulty lies in the extraction of the form factors from the three-point functions 
(P\V^\K) as a function of q 2 . One limiting factor is the cost of the calculation as one wishes to 
change the momentum transfer q, as generally this requires an additional propagator generation for 
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each q. Thus, to get both the shape and normalization of the form factors as a function of q 2 this 
can become an expensive calculation. 

One solution is to minimize the number of Dirac operator inversions that are required at the 
cost of introducing noise. Such stochastic methods are not new (for example, see Refs. j j40| , |4~I)]), 
and when implemented can easily save on the cost of the simulation. In particular, this was dis- 
cussed and preliminary results were shown for the shape of the form factor at this conference 
and this is an interesting approach to the problem. 

Another approach which can lead to more immediate phenomenological impact has been pre- 
sented by H. Na [ |42| ] for the HPQCD collaboration at this conference. Looking at the expression 
for the differential decay rate, we have 

where is the three-momentum of the outgoing kaon. What is then needed from the lattice is the 
particular form factor / + , which comes from the parametrization 



(K\sfc\D s ) 



Ui^ + ^-^^Hq 2 ). (5-2) 

q 



However, what is generally extracted from experimental calculations is the normalization of the 
form factors, or more specifically, |/ + (0)||V o |. Thus, a lattice determination of /+(0) is sufficient 
to determine \V CS \- 

This can be acheived with two simple ingredients. First, there is a kinematic constraint that 



states that the two form factors in Eq. ( pT2| ) are equal at q = 0: /+(0) = /o(0). Second, that the 



scalar form factor can be related to the vector form factor by 

qfl {K\sfc\D s ) = (K\sc\D s ) , (5.3) 
which leads to the important relation 

777 — 777 

MO) =/o(0) =-£ ^((K\sc\D s )\ j2=0 , (5.4) 

and the right hand side can be extracted with less noise than with the more direct method. Prelim- 
inary results, tested by looking at D s —* rj/v decays, were presented at this conference [fffil, and 
show a promising method to calculate the form factor at q 2 = with few percent precision. This is 
an essential goal for any method to calculate |V„|, or any other CKM matrix element. 



6. Nonleptonic decays 

When looking at decays of heavy-flavored states on the lattice, the primary focus is that of 
leptonic or semileptonic decays. The reason for this is that in these types of decays, there is either 
one or zero hadrons in the final state. One reason for this is simplicity: As you add hadrons in 
the final state, the calculation becomes more complicated and noisy. As such, fjj (zero hadrons in 
the final state) is much more simple to calculate than the form factors for D — > %iv, which has a 
single hadron in the final state. Of course, there is another reason for not looking at processes with 
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more than one hadron in the final state, and that is because it is not directly accessible. There is the 
Maiani-Testa "no-go" theorem that states that physical Minkowskian amplitudes are not obtainable 
from Euclidean correlation functions in an infinite volume p3|]. 

However, this is an obvious limitation in the field, as there is a wealth of information hidden 
within nonleptonic decays of hadrons. A large amount of effort, for example, has been spent trying 
to extract the amplitudes for K — » 2% decays (see for example Refs. Q ^5), |6[ 0, |||, |9|, 
even though it is limited by the no-go theorem. These various references apply various tricks and 
other techniques to bypass the Maiani-Testa theorem in order to attempt a lattice calculation of 

In fact, the decay of a kaon to two pions is complicated by several other features besides 
the Maiani-Testa theorem. It can occur in two channels, the M = 3/2 channel and the AI = 1/2 
channel. The former is more straightforward, but the latter is complicated by the addition of dis- 
connected diagrams and lower-dimensional operator mixing (leading to power divergences as the 
lattice spacing goes to zero), among other things. However, it is such an important quantity that 
there is an industry trying to circumvent these difficulties, so as to make any progress on a lattice 
determination of the K — > 2% amplitudes. 

In the same vein, there is a lot to be gained from studying nonleptonic decays of heavy-light 
mesons. For concreteness, I will focus on B — > Dn or B — > DK decays, and what can be gained 
from understanding these amplitudes The reason for this is primarily one of simplicity: These 
amplitudes, as I will discuss shortly, are free from several difficulties found in other processes, and 
they give essential insight into the CKM unitarity triangle. 

The effective Hamiltonian which governs these decays is given by [ ]52| ] 

H ^ = % I I [V: b V ui Cj(LL)Q'p cJ + V; b V ci Cj{il)Q ^ + h.c] . (6. 1) 

V 2 j=i t 2i=d,s 

The operators are denoted with j = 1 , 2 for the color mixed and color unmixed four-quark operators. 
What is most interesting is the ratio between the amplitudes for B — ► DK and B — ► DK (or similarly 
with K -> k). The ratio [0] 



Jt(B- 



0.1-0.2 



(6.2) 



allows us to extract the CKM angle y [|53|, |54|], if the various matrix elements [from four of the 
operators in Eq. ( |6.1[ )] can be determined from lattice simulations. 

These matrix elements are actually much simpler to calculate than those in K —> 2n decays. For 
one, there are no penguin contributions or disconnected pieces. In fact, of the four possible K —>2ti 
diagrams (in the AI = 1/2 case), only one enters in this analogous B — > DK decay. However, this 
does not remove the restriction imposed upon us by the Maiani-Testa theorem, and thus we must 
attempt to circumvent this problem. 

An initial, if rather crude, approximation is to look to Chiral Perturbation Theory (XPT), as 
has been done in the kaon case [55], and Heavy Quark Effective Theory (HQET) [56]. If one works 
out the heavy-light XVI, it is easy to show that at leading order in both m % and 1 /m^ c : 



(/sko^PT 



Jt(B s 



■if) 



Ji{B s -> D°) 



(6.3) 
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In this relation, I have used the fact that at leading order in HLZPT, both the b and c quarks are 



static, and thus the "decays" shown in Eq. < p.3\ ) are allowed via the weak operators in Eq. < \6.l\ ). A 
similar (and theoretically more accurate) relation can be found for the B — ► D% matrix elements. 

While valid at leading order, there could easily be large corrections coming from higher order 
effects. The corrections coming from HQET, where and m c are not infinite (but still degenerate) 
are most likely mild, as these corrections often can be absorbed into redefinitions of the various 
couplings [57]. The two corrections that are not likely to be mild are those coming from mt, ^ m c , 
and those coming from higher orders in the pion/kaon mass. 

For the first correction, where the heavy quarks are not degenerate, one can see immediately 
where this will be a problem. Setting the heavy quark masses to their physical values, the outgoing 
kaon or pion will carry a large momentum, well above the regime where XPT is reliable. In this 
case, one would not expect the determination of r% in Eq. (5.3) to be trustworthy. However, in the 
regime where my ~ m c , this is a reasonable approximation. 

As for higher orders in the light sector, often collections involving one-loop diagrams with 
kaons can approach the 15-30% level, depending on the quantity. While this is an unacceptable 
uncertainty for determinations of CKM matrix elements, one could hope that much of this uncer- 
tainty is cancelled in the ratio that is desired. Without doing the explicit calculation (both of the 
ZPT and the lattice calculation), nothing can be ascertained about this uncertainty, but it must be 
kept in mind before making quantitative claims. 

However, I would like to point out that this should not be taken as the ultimate approach to 
extract y and be able to make definitive claims with regards to the Standard Model. This is meant 
as a stalling point for this particular set of quantities. In a particular limit (nib = m c = °°,mu,d = 
m s = 0), this relation is exact, and allows for an initial approach to this problem from a lattice 
perspective. 

Ultimately, of course, methods to tackle the four-point functions themselves for B — ► DK 
and B — ► D% decays must be developed. These could include, but are not limited to, different 
(unphysical) kinematics to bypass the Maiani-Testa theorem or varying the boundary conditions 
(both of which have been used in the kaon sector). In addition to this, one must be able to get a 
handle on the final state interactions which are likely to be an important ingredient in these decays. 



7. Conclusions 



While in recent years much progress has been made in heavy flavor physics on the lattice, there 
is still quite a bit that needs to be done. As in the light quark sector, there needs to be numerous 
cross checks between different calculations using different techniques to ensure credibility in lattice 
calculations as a whole. At the same time, there have been steps made in new methods to better 
constrain quantities that have already been determined. Finally, I have proposed initial steps in 
lattice calculations of nonleptonic decays of heavy-light states, which are essential for a complete 
lattice picture of the heavy quark sector. 
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